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1. Introduction 

To name a date of begin of the ZFITTER project is difficult. The first papers on electroweak 
loop calculations by D. Baidin and O. Fedorenko date back to 1976, but in another context. In 
September 1983 the Dubna-Zeuthen group started activity, due to the begin of the four-year long 
stay of S. Riemann and T. Riemann at JINR, Dubna. The name ZFITTER was invented in 1989 
and replaced the former name ZBIZON for the software project. Finally, we chose the year 1985, 
when the article "Hunting the hidden standai^d Higgs" was published [1]. With this study, we began 
to take into account a finite, non-zero top quark mass in the radiative corrections, in the context 
of e+g^ -annihilation. To our knowledge, the paper contains the first plot confronting two LEP 
observables - weak mixing angle sin^ and Z boson mass Mz - with their dependence on the 
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unknown top quark mass and the also unknown Higgs boson mass Mh in the Standard Model 
[2-5]. We reproduce the plot here in figure 1.2, left. Both top quark and Higgs boson were not yet 
discovered at that time, and the actual experimental values for Mz and sin^ had too huge eiTors 
to be included into the plot [6]: Mz = 92.9 ± 16 GeV and sin^ dw = 0.23 ± 0.015. The numbers in 
the figure are based on the one-loop Standard Model prediction for Ar, the weak correction to G^, 
deserving few Unes of Fortran code. We remark as a curiosity that from 1985 to 2011, the article 
was quoted only once (by authors outside our group). 

The LEP/SLC collaborations made exciting measurements of the Z boson resonance and of its 
mass, width, weak mixing angle etc., with an unexpected final accuracy: 



Mz = 91.1876± 0.0021 GeV, (1.1) 

Tz = 2.4952 ±0.0023 GeV, (1.2) 

sin^^weak = 0.22296 ±0.00028, (1.3) 

sin^efftj = 0.23146 ±0.00012, (1.4) 

sin^el^s = 0.231 16 ±0.00012, (1.5) 

A^v = 2.989 ±0.007. (1.6) 



For the mass, this implies AMz/Mz ~ 10^^. For the various definitions of the weak mixing angle 
see section 10 of [8]. And Ny is the number of light neutrinos. 

Figure 1.3 shows the rise of accuracy for Mz due to LEP. Since the begin of the nineteen- 
nineties, a true scientific standard is the so-called blue-band plot of the LEPEWWG, based on 
ZFITTER [9-13] and another Standard Model package TOPAZO [14-16].' The March 2012 




Figure 1.1: ZFITTER authors: Lida Kalinovskaya, Pena Christova, Dima Bardin, Tord Riemann, 
Sabine Riemann, Andrej Arbuzov. Photograph made on occasion of 5th Helmholtz International 
Summer School "Calculations for Modern and Future Colliders ", July 23 - August 2, 2012, Dubna, 
Russia, CALC 2012. Right: Sasha Olshevsky, Arif Akhundov, Mark Jack. Fotos: © 2072 tordrie- 
mann@ googlemail.com; JINK, Dubna; A. Akhundov (priv.); Physics at FAMU. 

'Reference [11] (1992) appeared as a CERN preprint because, at that time, we considered this to be more presti- 
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Fig 1 Graph of sin^e^ versus Af^i "ifl^sriced by through 

radiative corrections The thickness corresponds to the range [^®^] 
30 GeV < < 40 GeV, the error bars indicate the accuracy 
expected at Z boson factories. 



Figure 1.2: Left: The first ever plotted LEP observables' dependence on the Higgs mass in 
the Standard Model ( reprinted from Physics Letters, A. Akhundov, D. Bardin, and T. Riemann, 
"Hunting the hidden standard Higgs", volume B166, p. Ill, Copyright (1986) [1], with permission 
from Elsevier. ) Right: Blue-band plot of the LEPEWWG [7] with a Standard Model Higgs boson 
mass prediction based on combined world data from precision electroweak measurements. 

version is reproduced in figure 1 .2, right. Both ZFITTER and TOPAZO are huge software packages 
with tens of thousands lines of Fortran code aiming at covering the complete known radiative 
corrections to the Z resonance peak in the reaction e^e^ — > //. The top quark was predicted by 
M. Kobayashi and T. Maskawa in 1973 [17] and discovered in 1995 with a mass of about 173 GeV 
[18, 19]. Top quark mass data from precision electroweak measurements and from direct searches 
are collected in figure 1.3, right. After discovery of the top quark, the LEP data were no more 
competitive. The agreement of direct measurements (in 'all data') and indirect measurements (in 
'all Z pole data') supports the validity of the Standard Model at the quantum loop level. Over the 
years the predictive power of the indirect searches for the Higgs boson mass improved considerably, 
and the discovery of the top quark was a crucial improvement for this. This is described in figure 
1.4. In 2012, the LHC collaborations reported the discovery of a scalar particle with a mass of 
about 125 GeV [20, 21], which fits into these expectations from the indirect searches. It might well 
be that it is the particle predicted by Peter Higgs in 1964 [22, 23]. 

We live with the ZFITTER project for about 30 years now, and ZFITTER is yet in use for a di- 
verse variety of applications, ranging from the global analyses of the LEPEWWG to many gradua- 
tion papers like e.g. [24]. Thirty years are a long term. It takes similai^ly long to prepare final results 
of big experiments at accelerators as LEP 1 (13 August 1989 - 1995), LEP 2 (1996 - 2 November 
2000), HERA (1992-2007). The final analysis of the LEP 1 data for two-fermion production was 
published in 2005 [25] by the LEP collaborations and the LEPEWWG, using ZFITTER v.6.42. 
The coiTcsponding enterprise for LEP 2 data was finalized these days [26], using ZFITTER v.6.43. 

The big laboratories invented scientific programs for a dedicated long-term preservation of 

geous than e.g. a paper in the journal "Computer Physics Communications" devoted to publication of software. It was 
submitted to the internet archive hep-ph in 1994. 
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Figure 1.3: Left: Z boson mass measurements at LEP. Earlier measurements are from UAl, UA2 
at SPS (CERN) (see text, not shown in plot) and from MARKII at SLC (SLAC). Right: Top quark 
mass measurements. 



the experimental data, under the label "ICFA Study Group on Data Preservation and Long Term 
Analysis in High Energy Physics" [27]. One might assume that this is a self-evident issue of 
any physics collaboration. Physics is the science of reproducible observations in Nature and of 
their explanations/descriptions, and reproducibility deserves storage. But long-term storage is an 
unsolved problem, worth of any (reasonable) effort. DESY, as an example, founded in 2009 a 
"DESY Data Preservation Project", mainly focusing on the HERA experiments HI, Hermes, ZEUS 
[28]. If such effort is justified for data, then it is also needed for the analysis tools, which were used 
for an extraction of the Model with its few parameters from the raw, or not-so-raw. Data. To our 
knowledge, the Big Labs do not plan to support long-term maintenance of software like ZFITTER. 
We, as authors, theoreticians and phenomenologists, have to mind by ourselves about maintenance 
of theory/phenomenology software. 

Everybody knows that the very details of a data analysis cannot be described by few words. 
But for precision studies they are truly essential. Sometimes we say: "The description of the 
program is the program itself." This is a helpful statement if "the program itself" is preserved over 
a long term in its state of use. ZFITTER did and does a lot to fulfil such a demand. See the webpage 
http : / /zf itter . com. 

Preservation demands effort. There are 17 people involved in the DESY Data Preservation 
Project. At the other hand, if a theoretician says: I care about the availability of my old software, 
people start to smile. This aim does not give true credit points for a scientific carrier, in what phase 
of the carrier ever. 

In fact, not only the so-called main author of ZFITTER, Dima Bardin, our "primus inter pares", 
tends to lose interest in active support of ZFITTER over the decades. This applies to all of us, 
mainly because of our interest in studying or inventing something new. Nevertheless, we collected 
in 2005 some volunteers into a ZFITTER support group, which submitted in that year ZFITTER 
v.6.42 and in 2008 ZFITTER v.6.43 [12, 13]. ZFITTER v.6.44beta dates in 2013 [29]. 
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Encouraged by the decreasing visibility of our ZFITTER support, in 2006 some experimen- 
talists tried to re-program in C-i~i- in a year's time the Standard Model library of ZFITTER from 
the published literature. Not just for fun, but in order to do better than ZFITTER: use a more 
modern programming language than Fortran, with more modularity than ZFITTER, a bit updated, 
with a GUI. In order to retain ZFITTER for a longer term. The project was proprietary until Au- 
gust 2012, and it faced two major problems. It proved to be impossible to do so without using the 
ZFITTER software itself to a large extent. Further, without cooperation with ZFITTER authors and 
the community of theoreticians, including extensive numerical cross-checks, such a project cannot 
succeed.^ 

Finally, there is much influence by institutes' directors and by the editors and publishers of 
physics journals on the engagement of scientists in the development of software. Not all of them 
seem to mind about proper acknowledgement and quotation of software. Some even say that soft- 
ware has no genuine scientific value by itself and advocate an absolutely free use of any software 
as common habit. If this would become common habit, nobody with inspiration and ambition 
would invest time to write complicated software for the use by other people, like the ZFITTER 
group - and other groups as well - does. We live in an academic world and we are valued by our 
scientific results, their originality, importance, curiosity, usefulness etc. Financing of our projects, 
of our working positions, our academic prestige depend on all that. We need proper quotation of 
our scientific results in case they are used. And we can only appeal and hope that the community 
understands this as a justified expectation, also for software. 

As a key feature of user-friendly support, we stored for many years all the relevant versions 
of ZFITTER at a webpage for anonymous download. We collected about three dozen versions, 
covering more than 20 years. There are coUegues who take the freedom to use ZFITTER as if it 
were open-source software in the strictest meaning of the word. Despite the facts that academic 
research deserves strict, proper quotation, and that there are licence regulations (for ZFITTER this 
includes the CPC licence). In some countries there are even legal regulations.^ 

It is the aim of these notes to give an overview on the ZFITTER project. Maybe they can help 
to see theoretical software in particle physics as an intellectual enterprise like the other inventions 
of physics research - experimental set-ups, data, hypotheses, models, theories. 

We would Uke to finish the introduction with two quotes. 

Several times we all thought that the ZFITTER project is in its final phase of dying out. See 
for example the remark of Dima Bardin at the symposium "50 Years of Electroweak Physics: a 
symposium in honour of Professor Alberto Sirlin's 70th Birthday", in the year 2000 [30]: 

"We would like to see the end of the ZFITTER project in the year 2000 and, therefore, a very natural question 
arises: What's next?" 

In the same year, members of the ZFITTER group were granted the prestigeous "JINR Award 
in Theoretical Physics" of the Joint Institute for Nuclear Research, Dubna, Russia. For a document, 
see here: certificate. The referee was Academician Prof. L. B. Okun from ITEP Moscow; he 

^See subsection 6.2. 

^Due to controversial positions, we closed the links for anonymous download from ZFITTER webpages in 20 11; in 
2012 the copies in the Andrew file system at CERN were removed. 
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Figure 1 .4: Higgs boson mass measurements. The upper limits and the fit values for Mh derive 
from a combination of virtual corrections to LEP and similar data, top and W mass measurements, 
performed by the LEPEWWG. The lower mass limit is due to LEP direct searches. The lower limits 
from data combinations are not shown. 

finished iiis estimate with the statement:'* 

"Overall, the project "ZFITTER Fortran program" represents a unique theoretical tool of world class. The project 
formed the basis of a close cooperation of experimetalists and theoreticians (with a series of workshops at CERN). With 
the accumulation of experimental data, the accuracy of the programs has been increased. The project has always found 
great interest at conferences. Its importance and the interest to it shows with numerous references in articles, reviews 
and monographies. 

In the long term, with the advent of more precise experiments, ZFITTER will allow to take into account all two-loop 
electroweak corrections. 

The series of theoretical articles on precision tests of the Standard Model at electron-positron colliders certainly deserves 
the award of the JINR prize 2000. 
Academician L.B. Okun" 

Our figures illustrate the development of mass predictions for Z boson (figure 1.3, left), top 
quark (figure 1.3, right), Higgs boson (figure 1.4). Here, ZFITTER has been useful until now. 
Okun's proposition that ZFITTER will be used also in future is being fulfilled. We can only hope 
that our write-up might help to convince the present particle physics community that ZFITTER is 
worth some support by now and in future. 

At the end of the introduction, we would like to reproduce the long(est) authors list of ZFIT- 
TER, see also http : //zf itter . com: 

A. Akhundov, A. Arbuzov, M. Awramik, D. Bardin, M. Bilenky, A. Chizhov, R Christova, M. 
Czakon, O. Fedorenko (1951-1994), A. Freitas, M. Griinewald, M. Jack, L. Kalinovskaya, A. 01- 
shevsky, S. Riemann, T. Riemann, M. Sachwitz, A. Sazonov, Yu. Sedykh, I. Sheer, L. Vertogradov, 
H. Vogt. Few of us are seen in figure 1.1. 

^The original document is in Russian, see statement. 
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The list is not complete. According to the conventions of the software library of "Computer 
Physics Communications", we should also include here all the co-authors who helped to prepare 
the program descriptions in 1989, 1999, 2005 [10, 12, 13]. 

2. ZFITTER in a nutshell, or: Is there a ZFITTER approach? 

We never used the label "ZFITTER approach". The reason is simple: There is no ZFITTER 
approach. If any, there is a kind of Dubna approach, or of Bardin's group's approach. 

Nevertheless, other people use this phrase. Let us collect some distinguishing moments which 
might be the origin of some popularity of ZFITTER, but also of one or the other of our scientific 
projects: 

• Unitary gauge. 

We are working in the unitary gauge when studying the renormalization of the Standard 
Model. Most of the other groups use the 't Hooft-Feynman gauge. But when looking at 
observable quantities, there is no difference left, due to the gauge invariance of perturbation 
theory. So, if everything is correct, there is no difference for the users. 

• On-mass-shell renormalization scheme. 

We are applying the on-mass-shell renormalization scheme, with few modifications. Other 
groups do the same for electroweak corrections. 

• Analytical treatment of QED corrections. 

ZFITTER is not a Monte-Carlo program. The Dubna group has an enormous experience in 
the analytical treatment of QED corrections, allowing us, sometimes, to come relatively close 
to the experimental set-ups by dedicated analytical integrations. Several different approaches 
may be chosen by users. The necessary computational time for fits to data is small compared 
to that of other projects. 

• Realistic observables and pseudo-observables. 

There is a plethora of observables, of quite different polarized and non-polarized cross- 
section combinations and asymmetries. Both so-called realistic observables (including real 
corrections) and pseudo-observables (after unfolding the realistic observables) may be used. 
With the different interfaces one may optimize a study appropriately. 

• Form factors. Modularity. 

We describe the effective Born cross-section in the Standard Model approach by (essentially) 
four (compex-valued) gauge invariant form factors per production channel.^ Plus a separated 
running QED coupling. This allows a modular programming, the efficient introduction of 
New Physics into the package, or the convenient export of the Standard Model corrections 
into another approach to the real corrections. 

• Higher-order corrections. 

Originally, we calculated the complete electroweak one-loop corrections to the Z resonance 

^Massive top quark production deserves six form factors [31 -33]. See also the in-depth discussion in [34]. 
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physics. By time, there became more and more electroweak, QCD, and mixed higher-order 
corrections available, and we had to implement them into ZFITTER. In the nineteen-nineties 
these implementations dominated our efforts for ZFITTER. It is not the genuin theoretical 
work we like, but has to be done. 

• Interfaces. Modularity. 

ZFITTER is not a fitting program. But from the very beginning, we were aware of the fact 
that a data analysis at e.g. LEP may rest on different sets of assumptions, being incompatible 
to each other. The notion of interfaces was developed. The interfaces call the kernel of 
ZFITTER with different compositions of input variables, real corrections, an effective Bom 
cross section. The users of ZFITTER can choose among few sample interfaces, or they write 
their own ones. 

• Flags. 

The use of ZFITTER may be controlled by flags to be set by users. Although this implies 
problems for the update by the authors, for users this is truly convenient. 

• Descriptions. 

ZFITTER is described for users at different levels of complexity. There are about 350 pages 
of instructions. 

• Simplicity of file structure. 

ZFITTER is easy to use. It has a simple file structure, is self-contained, and has a sample 
output. The installation at a computer is done and controlled within minutes. The installation 
of the user software, which is calling ZFITTER and performing data fits, writing tables and 
drawing figures, might be much more involved. 

• Numerical cross-checks. 

With very precise data available, as it was typical for LEP physics, a careful numerical con- 
trol of the theory software became mandatory. Here, a lot of coUegues, including competitors 
of ZFITTER, invested huge collective efforts. Without that, one could not trust the impres- 
sive physical results of that era, or the long-term reliability of the code. 

• Source-open programming. 

The scientific seriosity of ZFITTER is trustable because its source code is publicly available. 
Because we expect that the usual academic conditions of use are respected, notably the CPC 
licence, we say it is source-open software. The meaning of the word open-source software is 
controversial and it should not be used for ZFITTER. 

• Social aspects. 

A software package of some complexity, written for use by other people, must be supported 
and, in case, updated. The authors need some contact with the users. And, last but not least, 
some Ucence regulations have to be fixed if the authors want to get their academic credit, e.g. 
in form of proper citations. Since the authors of ZFITTER are employed at some institutions 
distributed over several countries, it is of vital importance that these institutions interfere in a 
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constructive way. We are happy that this did happen for a very long period, in view of several 
social restructurings of institutions and even countries. 

ZFITTER is a Fortran library of Standard Model predictions for the scattering process 

e^e-^ff (+7, +ny) (2.1) 

at energies in the range ^/s 20 GeV to 150 GeV, above quark bound states [meson factories] and 
below the top threshold. The package is to be called by interfaces 

• in the Standard Model; 

• in several model- independent approaches; 

• with Z' bosons and similar physics extensions; 

• etc. 

One may evaluate 

• realistic observables - polarized and non-polarized cross-sections and cross-section asym- 
metries with a variety of cuts on the final state 

• (pseudo-)observables like Mz, rz,a^°d, /?had, Ap^', A^, sin^e^^^, ••• 

• the form factors, for use in another analysis program, 
with different choices of input variables, e.g. 

• Mz, G^, m,, Mh, ttem, ••• 

• Mz, Mw, m,, Mh, a^m, "v, ••• 

3. Electroweak virtual corrections 

The first weak one-loop calculations were published as Dubna preprints by D. Bardin and his 
PhD student O. Fedorenko in 1978 [35-37]. Together with P. Christova, then also PhD student 
of D. Bardin, the by now famous articles on the complete on-mass-shell renormalization of the 
electroweak Standard Model were published in Nuclear Physics B [38, 39], for fermion scattering. 
See also reference [40]. The corresponding studies for weak boson production and fermion-boson 
scattering are unpublished [41, 42]. 

These calculations were complete, but assumed all fermions to be massless. When experiments 
showed that at least the top quark should be heavy, the top mass dependence was included [1, 43, 
44, 9]. Some studies of structural aspects in the renormalzation of the Standard Model are [45, 46]. 

All this was done in the unitary gauge, while the other groups usually worked with the 't Hooft- 
Feynman gauge. Later, this difference was of some value because an agreement of two calculations 
performed in truly quite different gauges establishes a powerful cross-check of the numerics. 

The first numerical program BFK (acronym for Bai^din/Fedorenko/Khiistova) was written in 
Fortran. 
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The Zeuthen partners, staying at Dubna from 1983 to 1987, worked out the renormalization 
of the electroweak Standard Model in the 't Hooft-Feynman gauge [47]. But because there was 
never a numerical program created, the results of this work were more or less useless; they had 
a mere educational aspect. Nevertheless, the experiences from that activity were used in order to 
perform the first calculation of flavor-changing Z boson decays into different lepton flavors.^ This 
was unpublished [48, 49]; see also [50]. An application to flavor- violating Z decays into different 
quark flavors was finally published [51]. Later, when we were working on precision predictions 
for LEP, the results could be easily transformed into the calculation of virtual top mass corrections 
in (flavor-diagonal) bb production at LEP and in Z decay [43]. And yet later, they were a starting 
point for studies of lepton number violation in e^e^ annihilation with heavy neutrinos [52] and 
with supersymmetry [53].^ 

3.1 Sirlin's approach 

The notion of form factors p and K in the weak neutral current were, to our knowledge, intro- 
duced by A. Sirlin:^ 

• p - contains the electroweak corrections to the Fermi constant 

• K- contains the electroweak corrections to the weak mixing angle sin^ By/ 

This approach allows to retain in the on-mass-shell renormalization scheme the Bom definitions 
also in higher orders: 



3.2 The HECTOR and ZFITTER approach 

For general 4-fermion scattering amplitudes, one needs a more general description. This was 
first introduced, to our knowledge, by the Dubna/Zeuthen group, in 1988, in the article "Elec- 
troweak Radiative Corrections To Deep Inelastic Scattering At HERA. Neutral Current Scattering" 
by D. Bardin, C. Burdik (Dubna), P. Khristova (Shoumen), T. Riemann (Zeuthen) [55]. The cor- 
responding software is retained until today as the Fortran package HECTOR [56]. So, strictly 
speaking, one might call this the HECTOR approach. 

*We mention for curiosity tliat tiie numerics of tliis one-loop project was performed with a pocket calculator TI-57 
with 50 program steps. The program had to be typed in after switching on. The price of the device was 120 DM in the 
CERN shop. 

^Several of the results in supersymmetry found in the literature turned out to be just wrong when we had a look at 
them. 

'^For a historical perspective, see reference [54]. 



sin-^ 0, 




(3.1) 
(3.2) 



where 




(3.3) 



(3.4) 
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We use four complex form factors p, Kini, Kfln, K"ini-fin for the parameterization of the weak 
amplitude, including the WW and ZZ box diagrams. In the article "A Realistic Approach to the 
Standard Z Peak" by D. Bardin, M. Bilenky, G. Mitselmakher (Dubna), T. Riemann, M. Sachwitz 
(Zeuthen) [9], we excluded the weak WW and ZZ box diagrams from the form factors, making them 
independent of the scattering angle. This is of advantage at LEP where these box diagrams have 
minor numerical influence. When form factors are independent of the scattering angle, analytical 
phase space integrations become possible. In ZFITTER, there is an option to switch between the 
approaches. 

The Bom amplitude is factorized into two pieces with vector coupling v, and axial vector 
coupling a, of a fermion / to the Z-boson: 

A(g)B= [uiy^{vi + aiYi)ui\ x [ufy^{vf + afY5)uf], (3.5) 

and is generaUzed by loop corrections to 

AvvY® Y+Aavyys «) y+AyaY® 775 +A«fl775 O 775, (3-6) 

or, equivalently, 

fiLL7(l + 75) «> 7(1 + 75) +ByLy® 7(1 + 75) +BLr7(l + 75) ® y+Byyy® 7. (3.7) 
With Z boson and photon exchanges: 

^ = J^y + ^z, (3.8) 
J^y ^ FA[y(^y\, (3.9) 
~ pz [775 ® 775 + Vqy® 775 + v,775 ® 7+ v^/7® 7] . (3.10) 

In Bom approximation, it is 

Vg/~v^xv/. (3.11) 
The form factors Fa, p, Kq, K;, jq,/ ai^e complex- valued functions of 5 and ?: 

= (3.12) 

= l + 5aQED{s), 

OCem = (3.13) 

af=l, f = q,l (3.14) 

Vf{s,tf^ = l-4sm^e„\Qf\Kf{s',t), f = q,l (3.15) 

Vql{s,tY^ = Vq+Vl-l + l6sm'^dw\QqQl\Kq,{s',t), (3.16) 

where we use 2e = — 1- From [57], eq. (3.3.1), we quote: 

= ieHlj;'hf^^p,f{s,t)!^y,{l + ys)'^y,{l + ys) (3.17) 
-MQe\4^e{s,t)y^(^y^{l + y5)-4\Qf\slKf{s,t)y^{l + y5)(S)y^i 
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The form factors may be used, in analogy to the Z decay matrix element of Sirlin, for definitions 
of effective vector and axial vector couplings and of an effective weak mixing angle: 

Gf = PefG^,, (3.18) 

sin^d^^ = sin^dw, (3.19) 

sin^d^^f = Kf sin^dw, (3.20) 

sin^e^^^^y = sin^dw- (3.21) 



The unique definition of an effective weak mixing angle is lost. 

The first applications of the calculations of weak corrections by the Dubna group were ap- 
plied, together with N. Shumeiko, to deep-inelastic scattering; see e.g. [58, 59]. There were close 
relations to the NA-4 experiment at CERN with JINR participation. The form factors p and K are 
simply related to the one-loop form factors introduced in the original renormalization articles by 
Bardin and Fedorenko (1978) [35-37] and Bardin, Christova, Fedorenko (1980) [38, 39]: 

Pef = l+F,,{s,t)-slAr, (3.22) 

= \+FQ,{s,t)-F,,{s,t), (3.23) 

Kf = l+F,^{s,t)-F,,{s,t), (3.24) 

K,f = l+FQQ{s,t)-F,,{s,t). (3.25) 

The corresponding relations of form factors Ftj and the Z boson matrix element are: 

X jr^ll + 75) «> 7m(1 + r5)F,, {s,t) - 4\Q,\sIy^ 8) 7m(1 + 75)^2. (^,0 

-4|G/|47m(1 + 75) ^ YnF^Q {s,t) + l6\QMf\strn ® 7M^pa (^'0 1 • 
So far we discussed matrix elements. The differential cross section for e^e^ — > // is: 

^ (l+cos^^?) [KT{r)+3ie{x{s) Kt{I)) + \x{s)\' MZ)] (3.27) 



J COST? 2s 

+ 2cosT^ [KFB{r) + ^e{x{s) Kfb{I)) + \z{^)\^ ^fb(Z)]|, 

with 

, , Of s 
X{s) = -4= 9 • (3.28) 

One has to caie about the choice of a constant Z boson width jz or an ^-dependent width Yz here 
[60]. 
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The effective couplings are: 

KT{y) = 

— Born 

Kt{I) = 

— Born 

Kt{Z) = 

— Born 

Kfb{7) = 
Kfb{I) = 

— Born 

Kfb{Z) = 

— Born 

Here, / denotes the initial state and / the final state. For the Drell-Yan process qq — )• /+/^, it is 
q = u,d and f = I. 

The Ccolor is the color factor, e.g. Ccoior = 3 for initial state quarks and final state leptons. 

A similar formula describes the special case of Bhabha scattering [61, 12, 62, 63]. The 1990 
article, a numerical comparison with W. Hollik [62], seems to be the most precise prediction for 
the effective Bom cross-section of Bhabha scattering until today. 

At the end of the subsection, we would like to emphasize that notions of form factors are 
not unique. We split, for purely phenomenological reasons, the matrix element into two pieces: 
a photon amplitude and a Z boson amplitude. The calculation of the running QED coupling is 
technically quite different from that of the weak loop diagrams. So this is reasonable. Gauge 
invariance justifies it, but only if handled with care. There are diagrams which mix a photon and 
a Z boson amplitude, and this is gauge dependent. So, in ZFITTER we decided to include all the 
corrections but the fermionic self-energy insertions, a bit arbitrary, into the Z boson amplitude. 

Such a separation of photonic and weak terms is wishful also for the charged current W boson 
mediated amplitude. But a gauge-invariant separation of (virtual and real) photonic corrections 
from W boson exchange is impossible. In HECTOR [64, 56], we found a way to do well-defined 
separations by considering logarithmic terms and just explicitely defining some rule. This really 
worked out. Years later, when building a softwai^e for e^e^ — )■ vvy, we could take over the weak 
charged current form factor into the Monte Carlo program of S. Jadach and Z. Was [65]. This 
reaction is, for v = Ve, unique: It depends both on neutral current and charged current amplitudes.^ 

Similar problems have been discussed when the ZFITTER form factors were adapted to atomic 
violation measurements [66]. 

3.3 Drell-Yan processes 

We went a bit into the details of a correct ansatz for the effective Born approximation in the 
Standard model in e+e^ -annihilation. The situation in a Drell-Yan process is quite similar. One 

^Bhabha scattering has also s- and f-channel exchanges, but only of neutral current type. 



Ccolor QjQf\Fy{s)\^ (3.29) 

2 2 
Ccolor Qi Qf, 

2Ccolor \QiQf\ Fy{s)*Pif{s,t)viVf (3.30) 
2Ccolor \QiQf\vB,iVB,f, 

Ccolor |P//(^,Ol'(l + + |V/|' + |V,-/|2) (3.31) 

Ccolor ivlj + ali)ivl f + alj), 

0, (3.32) 

2ccoior\QiQf\Fy{sypif{s,t) (3.33) 

2Ccolor \QiQf \ ClB,iaB,f-, 

2ccoior \pifis,t)\'^23i{viVf + Vif) (3.34) 

2Ccolor i2vBjaB.i) (2vB,faB,f) . 
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may study e.g. the running of the weak mixing angle sin^ B^^{s') as a function of the scale s' from 
a hard cross-section Oqi{s'): 



where both hard scattering cross-sections Oq{uu — )• l^l^) and ao{dd ^ l^i ) depend on four com- 
plex valued, process-dependent form factors Pqi, Kq, Ki,Kqi with q = u,d. The Gq depends on s', but 
also on the scattering angle 6. Further, we have not only initial and final state photonic corrections, 
but also initial-final state interferences. 

An elegant way to cover at least part of the complexity of all this in a modern QCD Monte 
Carlo program is the following: 

- Define a photon exchange amplitude. 

- Define a Z exchange amplitude. 

- Split the Vqi into a Z-part and a photon part: 



- Assume a Born like structure with form factors p,Vq,vi and put the deviation from that 
structure, which is contained in the difference [vqi —VqVi), into the photon amplitude. 

In an unpublished paper of 1991 [67], A. Leike and T. Riemann worked out the influence of 
Z'-physics on the evaluation of weak form factors. The idea of reshuffling matrix elements in form 
factors was invented there and is now independently re-used as clever inclusion of ZFITTER's 
weak form factors into a Monte Carlo code, which was originally made for the description of QCD 
connections to Drell-Yan processes.'*' 

Evidently, once there are accurate data, one has to carefully understand how to model the 
correct physics ansatz with a smaller number of parameters. This is under study by experimentalists 
presently. 

4. Real corrections, mostly due to QED 

Around 1983 we began to envisage some contribution to the description of the Z boson reso- 
nance as it was planned to be studied at LER There existed several articles on electroweak radiative 
corrections. Let us mention the electroweak study by Wetzel in 1982 [68] and that by Lynn and 
Stuart in 1984 [69], or the MC program MUSTRAAL by Berends, Kleiss, Jadach in 1982 [70]. 
It was not evident to us that we might contribute some novel results, and we decided therefore to 
study real photon emission first. 

The Dubna group has an enormous experience in the analytical treatment of QED corrections, 
first mostly applied to ^-channel exchange processes. We should mention here the close contacts 
with Dubna experimentahsts of the NA-4 collaboration at CERN. The subtraction method for the 
treatment of infra-red singularities was worked out in 1976 in a seminal paper [71]. The divergent 
part of the cross-section is, in simplified form, integrated over the whole phase space, and at the 
same time subtracted from the exact squared matrix element. The difference can be integrated 

"'W. Sakumoto, private information. 



ao(/)=^„ Oo{uu^l^l ) + Ooidd^l^l ) 



(3.35) 




(3.36) 
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numerically, and the isolated term is sufficiently simple for an analytical treatment. In practice, this 
can become quite involved, see reference [72]. 

The first articles treated just photonic corrections, taking into account mass effects. The very 
first one was on pure QED corrections in e^e^ annihilation, by A. Akhundov (Baku), D. Baidin 
(Dubna), O. Fedorenko (Petrozavodsk), T. Riemann (Dubna): "Some Integrals For Exact Calcula- 
tion Of QED Bremsstrahlung", an unpublished JINR Dubna preprint [73], followed by [74, 75]. 
Then we extended the integrational technology to experimental set-ups with Z boson resonance 
phenomena, including mixing phenomena of Z boson and photon. This sounds easy, but there were 
several conceptual problems to be solved. As a result, ZFITTER relies now on several versions 
of semi- analytical formulae with low-dimensional numerical phase space integrations left. At the 
time of LEP experiments, this was extremely useful. For an unfolding of measured cross-sections 
into pseudo-observables, or for multi-dimensional fits, the computing time of an analysis code was 
absolutely decisive. The inclusion of certain kinematical cuts was a wish expressed by experi- 
mentalists. Computers were not so advanced. There were no personal computers, and workstations 
were also not yet on the market. In Dubna, there were one or two terminal stations for theoreticians, 
and we had to queue up every day. In Russian Winter, the terminal room (with one terminal) was a 
bit cold at temperatures close to zero centigrades, because the windows did not close exactly. The 
upper left corner of the terminal screen was blind. Often the terminal in the theory building was 
blocked by Riemann, Bardin, Akhundov from 9 to 12 in the morning. Not everybody was amused. 

In case of quark-pair production, or Z or boson decays into quarks, the final state will 
get QCD modifications. The corrections are contained in so-called radiator functions. Their im- 
plementation in ZFITTER relies on calculations by a variety of coUegues and is described in the 
various ZFITTER descriptions, notably in references [9, 11-13]. Useful representations are also 
e.g. [76-79]. 

The treatment of the complete set of QED corrections related to real emission of photons in 
ZFITTER is quite specific. The higher-order corrections have been typically taken over from the 
Uterature, as it is documented, notably in references [12, 13]. An important example is reference 
[80]. The main work had to be performed at one-loop order, plus soft photon exponentiation. It 
was clear that the numerical effects will be important for the experimental analyses. There were 
several Monte-Carlo programs available, e.g. [81, 70, 82-85] and the references therein. See also 
the report [86]. We aimed at an alternative, analytical integration of the three-dimensional photon 
phase space integrals. The necessary techniques have been developed step by step over a longer 
period, and originate to a large extent from studies for deep-inelastic scattering, e.g. IN — )• IX [72]. 
In the presence of the Z boson resonance in the 5-channel, one is faced with the additional need to 
perform a correct treatment of the Breit-Wigner propagator, a truly complex function. Further, there 
is a mixing of photon and Z boson exchange. This yZ mixing was studied e.g. in [47, 51, 87, 88]; 
this issue was settled by a formal Dyson summation of the 7, Z propagator matrix. The Z boson 
propagator with the finite width may become an issue for analytical integrations. In squared matrix 
elements, we are faced with 77, yZ and ZZ interferences. The latter are dominating around the 
Z boson pole, and they will contain squared Z boson propagators. To perform analytical phase 
space integrations with such a term inside looks difficult. An important, simple idea is to perform 
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a partial fraction decomposition in order to linearize the integrand: 
1 

s-Ml + iMzTz 

MzTz ^\s-Ml + iMzTz 

At first glance this looks bizarre because the complete answer seems to carry an overall factor 
s/{MzTz)- Evidently, one may use complex integration theory, so this is good. The overall pre- 
factor gets divergent for vanishing Z width, but this is the technical expression of the well-known 
radiative tail, so this is also good. 

We tried the approach, and calculated the complete one-loop QED corrections for the total 
cross-section and the forward-backward asymmetry around the Z resonance without a cut. The 
results for initial state radiation, final state radiation and the initial final state interferences were 
rather compact and looked explicitly reasonably behaving. The results were published as preprint 
in [90] and refined a bit in [91]. The paper could not be published in Nuclear Physics B because 
the referee found it not close enough to the experimental set-up. Nevertheless, it is a nice piece of 
work and served for many years as an important numerical etalon for precision comparisons. 

As a by-product, we understood that one may calculate the photonic corrections to the initial- 
final state interference of the yZ interference as the arithmetical means of the corrections to the ZZ 
and 77 initial-final state interferences: 

Rini-fin{Z,Z2) = ^[Rini-finiZ,Z)+Rini-fin{Z2,Z2)], Z2 = Y- (4.2) 

Here, Z2 is a second vector boson. For a proof see reference [87] . This is not of utmost importance 
here. When we later studied QED corrections for Z,Z' production with a heavy Z' boson, then we 
had the newly appearing initial-final state part of the ZZ' interference at the disposal without a new 
calculation [92-96]. 

Later we refined the techniques, and finally ZFITTER enables the calculation of 



2/Mzrz Ks-M^ + iMzTz s-M^-iMzTz 



2^—TT^\ (4-1) 



• exact, completely integrated one-loop photonic corrections without cuts [91]; 

• convolution integrals for cross-sections with soft photon exponentiation [97] ; 

• the corresponding angular distributions [98]; 

• convolution integrals with integrated angular cuts [99]; 

• convolution integrals with integrated acceptance cuts, combined further with an acoUinearity 
cut [100]. 

"in fact, it took us nearly half a year of heavy fighting with SCHOONSCHIP, because we did not agree, at the Z 
boson peak, with the numerics of the Monte-Carlo program MUSTRAAL [70, 89]. The MUSTRAAL was available 
via CPC, and we could run it at Dubna. The mistake was, as often, trivial, but influential. The final 5 digits agreement 
convinced us that our Breit-Wigner treatment makes sense and is operational. 
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The sophisticated final state phase space treatment with cut on the acoUinearity final state fermions 
goes back to G. Passarino (1982) [101] and is relatively close to realistic experimental cuts for 
lepton final states. The complete analytical QED coiTcctions were worked out for this case by 
M. Bilenky and A. Sazonov [100] and became part of ZFITTER. The truly nice paper remained 
unpublished, unfortunately. Later, we recalculated these connections for ZFITTER from the scratch 
(unpublished, see references [102- 105]). We performed two minor corrections and got very nice, 
compact formulae for the special case of no cut for the fermion production angle [106]. 

Finally, all this was sufficiently close to what the experimentalists could derive from their 
Monte-Cai^lo simulations for a confrontation with theory. 

We wrote relatively monstrous programs in Veltman's SCHOONSCHIP [107, 108] to be run at 
a CDC-6500 main frame at JINR Dubna. Bardin and Fedorenko where, in parallel with Vladimirov 
and Tarasov, among the first using SCHOONSCHIP at JINR in 1976.^^ Collegues from Moscow 
came to JINR regularly in order to use the CDC-6500 main frame because comparable computers 
were subject of the US embargo policy and thus not available for civil use in Soviet Union at 
that time. JINR, Dubna, as an international research center, was privileged in that respect. A 
comprehensive review on the use of computer algebra at JINR is [109]. Later FORM [110, 111] 
was invented by Jos Vermaseren and we could run it at personal computers. The first article typeset 
in latex was presumably [98], and the first article submitted to the hep-ph archive dates in 1994; 
The archive hep-ph was opened in 1992. 

At a certain moment we realized that analytical integrations are fine; but if the sensitivity to the 
Z boson width becomes sufficiently large, then it will matter whether the width is a pure constant 
Yz as in a normal Breit-Wigner function, or whether it arises from a quantum field theoretical 
calculation and will thus depend on the kinematics, Tz{s). In the latter case, it is (roughly speaking) 
the imaginary part of the Z boson self-energy function, which is by itself ^-dependent; and for initial 
state corrections i^'-dependent. The s' is one of the integration variables. We remembered that the s- 
dependence is, to a very high accuracy, justMz '^z{s) = (s/Mz) Fz, and this observation enables us 
to change the propagators into functions with a constant width, allowing not only a good estimate 
of the different approaches, but also furtheron the analytical integrations: The differences of mass 
and width in the two approaches derives from the foUwoing identity [60] : 

^ TT = c 1-^ , (4-3) 

5 - M| + iMzTzis) s-m^ + imzYz 

with 

mz = Mz — = Mz - 34 GeV, (4.4) 

2Mz z , 

r3 

72 = Fz ^ = Fz - 0.934 MeV - 1 MeV. (4.5) 

^ 2M| 



'^A. Akhundov, D. Bardin, L. Bobyleva, V. Gerdt, I. Shidkova, W. Lassner, V. Rostovzev, O. Tarasov, R. Fedorova 
and D. Schirkov received in 1986 the JINR Award in Theoretical Physics for "Introduction, development and use of 
computer systems for analytical calculations at central computers of the central computing installations of JINR". We 
are grateful to V. Gerdt for a clarifying email exchange. 

'■'We are grateful to Andrei Kataev reminding about this fact. 

'^The Z boson mass shift was discovered by a numerical study of the Z boson peak in parallel in [1 12]. 
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Here, Mz = 91187.6 GeV and Tz = 2.4952 GeV have to be chosen as the usual PDG values. 
Later we worked out an approach to a model-independent Z boson peak analysis inspired by S- 
matrix theory, relying naturally on mz, Yz- Not only for the Z boson peak cross-section, but also for 
asymmetries. The point here again is a proper treatment of QED corrections [113-1 18].'^ 

In fact, the idea to use mz,yz instead of Mz,Tz was born while listening to a talk on string 
theory at a conference, while reading a paper on QED coiTections with complicated phase space 
cuts by Passarino [101]. 

The Z boson parameter relations (4.4) and (4.5) become essential when two-loop electroweak 
corrections are determined in ZFITTER. This is carefully described in [119], where the complete 
electroweak two-loop conections to the leptonic weak mixing angle have been calculated. See also 
section 6. It is remarkable that the shift of the Z boson width due to the change of scheme {s- 
dependent or constant Z boson width) amounts to 1 MeV and is larger than the conesponding shift 
from the genuine weak NNLO corrections. Compared to the experimental en^or of 2.3 MeV, the 
shift is small. The authors of [1 19] did not take the correction into account because it is formally 
beyond the NNLO order and thus among the systematically neglected terms. One should consider 
the term as an indication of the size of unknown higher order terms. 

What we describe here is about the state of real emission affairs in ZFITTER at the end of the 
nineteen-eighties. Final state mass effect treatments were refined in [120- 122]. Some additional 
QED corrections, due to light fermion pair emission and higher order photonic effects, needed for 
a proper treatment at LEP 2 energies were later added [123, 124]. See also reference [125]. 

Careful studies of ZFITTER physics updates originated in these years [126, 126, 78, 127]. 

5. Competition and cooperation 

5.1 1989 - First LEP publications 

In 1989, the world changed quite a bit. Participation at the Ringberg Workshop on LEP physics 
in Germany became possible [128]. The NATO supported RADCOR conference on radiative cor- 
rections and their applications to experiments in Brighton, the first one of a series, was open to 
Eastern Country physicists [129, 130]. We remember the stimulating athmosphere of the 1989 
LEP physics workshop at CERN [129, 131]. And LEP became operative in August 1989. The 
first months were exciting. A good knowledge of radiative connections was needed from the very 
beginning, just in order to discriminate between trivial radiative effects and New Physics. Several 
unpublished ZFITTER related theory studies appeared in this period, e.g. [132-135]. In [135], 
approximate parameterizations of 0{aas) connections [136] were derived in order to speed-up the 
numerics. The Fortran routines of B. Kniehl [137] improved this later furtheron. The LEP collab- 
orations performed the first Z line shape analyses. We were closely related to the L3 collaboration 
[138-145] and to DELPHI [146- 150]. A review of the latter is [151]. 

Among the first DELPHI papers was [146]. From the ZFITTER group, D. Bardin and G. 
Mitselmaker were DELPHI authors. The paper quotes for the theory on the Z line shape G. Burgers 
[152] and A. Bonelli et al. [153]. In [149], the Z line shape analysis used the softwai^e packages 

'^The corresponding software package SMATASY is supported by Martin Griinewald. 
'^Ayres Freitas, private information. 
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ZAPPH and ZHADRO by G. Burgers [152]. In [148], March 1990, our papers [9, 100] are quoted. 
And in [150] the package ZFITTER/ZBIZON with reference to the internal note DELPHI 89-71 
PHYS 52 and to [9, 10] was used." 

A similar approach was observed in the L3 collaboration, were ZFITTER authors T. Riemann, 
M. Sachwitz and H. Vogt were collaborating in 1989. The internal note L3-001 [138] quotes G. 
Burgers [152] and CERN 89-08, but also our paper [60]. The Z line shape analysis seems to be 
based on papers by Cahn [154] and Borrelli et al. [153]. In [139], internal note L3-003, our package 
ZBIZON is quoted with reference to L3 Internal Note 679 as well as [60] and the Zeuthen preprint 
PHE 89-19 [98]. Back-up radiative corrections had been studied with ZBIZON. For the very Z line 
shape fits they used again Borrelli et al. [153], Cahn [154], and a paper by Jadach et al. [83], for 
Bhabha scattering. In [155], internal note L3-004, the paper on the Z boson parameters [60] was 
quoted. 

A bit later it became more and more common to use ZFITTER in DELPHI and L3, but also 
in OPAL. While ALEPH used the package BHMAVOH by F Berends, M. Martinez, W. HoUik et 
al. [156, 76]. We mention these very first papers on LEP physics results because they demonstrate 
that there was a true competition of the analysis packages and our ZBIZON/ZFITTER package was 
accepted step by step, but not from the very beginning. 

5.2 1992-2012 - LEPEWWG and global fits 

The LEP Electroweak Working Group was founded in 1993.'^ Soon after the first measure- 
ments at LEP the quest was expressed for combined data analyses with a fourfold statistics com- 
pared to a single experiment. Originally a group with members of the four LEP experiments, led 
by Jack Steinberger, investigated the combination of the Z line shape [157]. In 1993 Dorothee 
Schaile was asked to take over the coordination of the group and she had then already ideas on the 
inclusion of other electi'oweak observables into a combined analysis. They called themselves the 
LEP EWWG. The first publicly accessible document with this name is also the initial summary of 
the LEP results for the electroweak Summer conferences in 1993, which then appeared annually 
[158 - 160]. The LEP EWWG was lead by D. Schaile from 1993-1996. When she became pro- 
fessor in Munich, Robert Clare took over the coordination of the LEP EEWG.'^ The present chair 
is Martin Griinewald. The final paper on LEP 1 data appeared in 2005 [25], nearly a decade after 
closing LEP 1 in 1996, while the analysis of LEP2 data (finalized data taking in 2000) was finished 
these days [26]. 

The ZFITTER group members, as well as the authors of other physics software packages used 
by the LEPEWWG are not members of the LEPEWWG. They are consulted in case. 

5.3 1995 - The Electroweak Working Group Report 

The work of the LEPEWWG and of the four LEP collaborations relied on ZFITTER and 
TOPAZO, and also on the BHMAVOH package, and on many other resources. Because of this 
role of establishing a kind of world standard, the community felt the need of cai'eful numerical 

'^ZBIZON is the former version of ZFITTER. 
"^We are grateful to Dorothee Schaile for private information, 
are grateful to J. Mnich for a clarification. 
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checks on their predictions. One is confronted with multi-parameter problems, different calcula- 
tional schemes, some freedom of input choices, in the presence of approximations and dedicated 
omissions, of misunderstandings and, sometimes, mistakes. 

At a certain moment, the community has to set benchmarks. The result of a year-long 
workshop is the collection "Reports of the working group on precision calculations for the Z 
resonance" , edited by D. Bardin, W. Hollik, G. Passarino. It was published as a CERN Yel- 
low Report, CERN 95-03 (31 March 1995), http://cdsweb.cern.ch/record/280836/ 
f iles/CERN-95-03 .pdf. 

Part of this document is the "Electroweak Working Group Report", which was two years later 
submitted to the arxive/hep-ph [76].^° This work is one of the basics for the successful work of 
the LEP Electroweak Working Group. It is until now one of the most important collections of 
Standard Model higher order corrections for e+e^ -annihilation. 

5.4 Higher order corrections in ZFITTER 

During the 1995 CERN workshop and shortly after, a lot of additional higher-order corrections 
were calculated and included into ZFITTER. We give here just a (presumably not complete) list of 
references and refer for any detail to the ZFITTER descriptions: [77, 88, 119, 161, 162, 137, 163 - 
167]. Later, further improvements were added [168-177]. 

Until now, we did not yet include into ZFITTER the existing parameterization of the rather 
small bosonic two-loop weak corrections to the weak mixing angle [175]. The fermionic correc- 
tions are covered, as well as the complete weak two-loop corrections to the W boson mass. For a 
complete treatment of the weak two-loop corrections to the Z boson width, the corrections to the 
form factor pz are lacking yet. For this reason, the quite good agreement of the higher-order ap- 
proximations to Tz with the so far known pieces of the complete two-loop result are an indication 
that the final answer will be close to what we have already. 

Generally speaking, we try to control about four to five digits of the predictions aiming at 
such a physical theory precision. One quote from the report [76] is interesting because it sheds 
some light on the progress of the so-called technical precision (precision under fixed, maybe not 
realistic conditions): " ... compare results of independent calculations. Such a comparison has 
been done once for Ar, and an agreement of up to 12 digits (computer precision) was found [14]." 
Ref. [14] was private communications of D. Bardin, B. Kniehl and R. Stuart in 1992. This has 
to be compared to a three digits agreement between two Bhabha cross section calculations in a 
comparison, performed few years earlier in 1990 [62]. Later, in 2002, a precision of up to 12 digits 
was reached in practice for complete virtual one-loop calculations, and of 5 digits with inclusion 
of real corrections [32, 178, 179]. 

6. ZFITTER 2013 

6.1 From ZFITTER v.6.42 to ZFITTER v.6.44beta 

The most recent publicly available ZFITTER version is ZFITTER v.6.43 (17 June 2008) 
[12, 13]. It agrees with ZFITTER v.6.42 up to a coiTcction of a non-influential typo and was 

^^Now it is also available as a pdf file at CERN, in CERN 95-03. 
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released by the ZFITTER support group (A. Arbuzov, M. Awramik, M. Czakon, A. Freitas, M. 
Griinewald, K. Monig, S. Riemann, T. Riemann, see http : //zf itter . com). The ZFITTER 
group was reorganized in February 2012 and consists now of A. Akhundov, A. Arbuzov, D. Bardin, 
R Christova, L. Kalinovskaya, A. Olshevksy, S. Riemann, T Riemann. 

Recently, we have included into ZFITTER v.6.44beta (20 January 2013) the final results for 
the <^{cc^) QCD corrections to the Z-boson and W-boson quarkonic partial widths and to the so- 
called /?-ratio by R Baikov et al. [29, 177]. As may be seen from figure 6.1 and from table 
6.1, the numerical shifts in the widths amount to less than 0.3 MeV and are thus well below the 
experimental errors, e.g. at LEP or at an anticipated GigaZ option of an ILC [180].^' A fit formula 
for the complete electroweak two-loop corrections to the W-boson mass [171] was already included 
in ZFITTER v.6.42. The final exact results for the complete electroweak two-loop corrections to 
sin^ for light fermions / [119] and the two-loop electroweak fermionic corrections to sin^ 6^^ 
[176] have to be included yet into ZFITTER. They are known to be small corrections compared 
to the fit formula [174] covered in ZFITTER since v.6.42. Already these corrections are small 
compared to the present experimental errors for the gauge boson widths, see table 6.1. For the 
leptonic weak mixing angle, they are of the order of the experimental en^or: Compare the Particle 
Data Group value of (1.4) with the last row in table 6.1. The comparison shows even a systematic 
deviation of the two values. This deviation traces back to the handling of the hadronic contributions 
to the photonic vaccum polarization. Changing the ZFITTER default by flag setting ALEM=2 
into a variable input and setting this to Aa^J^{Mz) = 0.02750, produces a shift of the ZFITTER 
prediction towards the PDG value.^^ See the changes shown in table 6.2. Just to mention, the 
influence of Aal^J^{Mz) on the Higgs mass prediction is visualized in figure 1.2, right. Here it is of 
minor importance, but visible. 

Presently, there are controversial positions concerning ZFITTER's 'conditions of use' and 
the ZFITTER software licence http : //cpc . cs . qub . ac . uk/licence/licence . html 
granted to the authors by Elsevier's Computer Physics Communications Program Library - Pro- 
grams in Physics & Physical Chemistry. For some details see http://zfitter. com. Until the 
issue is settled, actualized versions of ZFITTER will stay at the beta level and cannot be released. 

Sooner or later, the LHC is becoming a precision tool and the community feels some steady 
need of high-precision Standard Model predictions. Both for use in global fits and for specific 
cross-section predictions, notably of Drell-Yan processes via the Z resonance. This need would 
become even more pronounced if the ILC project would substantialize [180]. 

Regrettably, we see today no alternative project to ZFITTER in the field of precision Standard 
Model predictions. In the mid-nineteen nineties there were three competing (and cooperating) 
projects at the disposal [76]: BHMAVOH by W. Hollrk et al., TOPAZO by G. Passaiino et al., and 
ZFITTER by D. Bardin et al. BHMAVOH was available on request, and the latter two aie publicly 
available. To our knowledge, updating and user support have been minimized for TOPAZO and 
BHMAVOH [156]. 

^' A detailed numerical study is in preparation [181]. 

^^Taking into account the uncertainty Aa^J^{Mz) = 0.02750 ±0.00035 [182], the corresponding predictions in table 
6.2 vary: rz(/l+/i") by ±6.7 X 10"^ , Tz by ± 1.2 x lO""* , rw(/v) by ±2.2 x 10""* , r^/ by ± 2.2 x 10""^, Mw 
by ± 7.5 X 10^^ , sin^ 0^^?' by ± 5.0 x 10^^. The latter is about the value of the the experimental error. 
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Figure 6.1: The influence of the i&{ccj) QCD corrections [177] on the W and Z boson widths. 



6.2 A comment on the Gfitter project 

Sometimes the Gfitter project is considered as an independent implementation of Standard 
Model predictions for some pseudo-observables, and as a true scientific alternative to ZFITTER (for 
these pseudo-observables). We do not share this opinion and would Uke to give a short, clarifying 
comment on the situation. 

The Gfitter project was started in Summer 2006 and presented to the public in De- 
cember 2007, at the kick-off meeting of the German "Helmholtz Alliance for Physics at 
the Terascale", see the slides at http://indico.desy.de/materialDisplay.py? 
contribId=36&sessionId=15&materialId=l&conf Id=477. Until August 2012, the 
Gfitter software was proprietary, but by private information^^ it became known that the Standard 
Model library of Gfitter, Gfitter/GSM, was relying on the FORTRAN package ZFITTER v.6.42 
and was created to a large extent by copy-paste-adapt. Without any proper citation in the academic 
meaning of the word. 

There are several versions of the program Gfitter. 

• Gfitter/GSM (Summer 2006 - July 2011) relies essentially and directly on the Standard 



■^^Private information from and documentation by A. Akliundov, S. Riemann, T. Riemann, March to May 2011, 
http: // z fitter . com. Further, a German ombuds person's report announces in July 2012 : "A diploma thesis 
derives from ZFITTER in the sense that 8200 lines have been taken over by copying from ZFITTER." In the thesis work 
the kernel of the Gfitter/GSM software was written (in collaboration with others), and its text delivered basic building 
blocks for the so-called main article on Gfitter [183]. A third evidence for the confidential take-overs may be found 
in the unpublished version of Gfitter of July 2011, where about 100 to 200 identities are denoted, by the Gfitter/GSM 
authors, to originate from ZFITTER v.6.42. On occasion of the Erratum [184] to [183], ZFITTER authors wrote a letter 
to the Editorial Board of "European Physical Journal C" (14 September 2012); see letter-to-the-epjc-editors.pdf. 
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Table 6.1: ZFITTER 6.44beta, with the input values a, = 0.1184, Mz = 91.1876, Mh = 125, m, = 173. 
The dependence on elect roweak NNLO corrections is studied for IM0MS=1 (input values are aem,Mz,Gn). 
AMT4= 4: with two-loop sub-leading corrections and re-summation recipe of [23-28] of [13]; AMT4=5: 
with fermionic two-loop corrections to Mw according to [29,30,32] of [13]; AMT4=6: with complete two- 
loop corrections to Mw [37] and fermionic two-loop corrections to sin^ 0^'"''''^ [52] of [13]. IBAIKOV=0 
(no QCD corrections) or IBAIKOV=2012 [177]. 



AMT4 


4 


5 


6 


Diff. 


Exp. Err. 


1BA1KOV=0 




83.9782 


83.9748 


83.9807 


0.0059 


0.086 




2494.7863 


2494.6019 


2494.8688 


0.2669 


2.3 


rH^(/v) 


226.3185 


226.2877 


226.2922 


0.0308 


1.9 




2090.3308 


2090.0465 


2090.0882 


0.2843 


42 


Mw 


80.3578 


80.3541 


80.3546 


0.0037 


0.015 


sin^ el'^^ 


0.231722 


0.231791 


0.231670 


0.000121 


0.00012 


1BA1KOV=2012 




83.9782 


83.9748 


83.9807 


0.0059 


0.086 


Tz 


2494.5591 


2494.3747 


2494.6416 


0.2669 


2.3 




226.3185 


226.2877 


226.2922 


0.030 


1.9 


Fw 


2090.1117 


2089.8274 


2089.8691 


0.2843 


42 


Mw 


80.3578 


80.3541 


80.3546 


0.0037 


0.015 


Sin 


0.231722 


0.231791 


0.231670 


0.000121 


0.00012 



Model implementation of the ZFITTER softwai^e. On top of that, Gfitter/GSM contains 
few add-ons. The electroweak add-on of Gfitter/GSM, compared to ZFITTER v.6.42, are the 
bosonic two-loop connections to the weak mixing angle in Amwramik et al. [119]. They are 
small; see the discussion above. The complete two-loop parameterizations in [119], in turn, 
have been made with use of ZFITTER v.6.42. As a consequence, it is formally correct to 
quote for the parameterization only [119], but one should have in mind that there is inside 
also ZFITTER. There is also a QCD add-on of Gfitter/GSM (201 1), compared to ZFITTER 
v.6.42 (2006), based on [185]. It is also numerically small (see the discussion above) and is 
implemented in ZFITTER v.6.44beta. 

Use of this Gfitter version deserves a citation not only of [183], but also of [12, 13], for using 
ZFITTER v.6.42, according to its licence. 

• Gfitter/GSM (August 2011 till August 2012) relies on a proprietary implementation of Stan- 
dard Model corrections which were based on a parameterization tracing back to Cho et al. 
(1999) [186], which in turn is based on an electroweak one-loop calculation published in 
1994 [187]. There have been made improvements later, and in a recent article by Cho et al. 
(2011) [188] the authors confirm the reliability of their pai^ametrization by comparing them 
with ZFITTER v.6.42 predictions. These parameterizations ai^e used in Gfitter furtheron, and 
overlayed with the most recent higher-order corrections mentioned. 
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• Gfitter_1.0 has been released publicly in September 2012. The Standard Model library Gfit- 
ter_1.0/gew relies presumably on the same parameterizations as Gfitter/GSM (2011). 

The different versions of Gfitter rely in one way or the other on ZFITTER v.6.42. We further 
remark that without studying the numerical reliability of Gfitter, to four or five significant digits, 
the scientific value of the inclusion of NNLO weak and QCD corrections in Gfitter remains 
questionable. According to our standards, Gfitter simulates Standard Model predictions with un- 
known precision. It is a nice tool for the production of figures for the illustration of Standard Model 
physics. Possibly it is useful for studies beyond the Standard Model. 

7. Conclusions 

A talk on history and features of the ZFITTER project was presented at LL2012, the eleventh 
"Loops and Legs" meeting. Its title was "ZFITTER - 20 years after". The "Loops and Legs" 
conference was founded by the Zeuthen Theory Group in 1992 when the Zeuthen Institute for 
High Energy Physics of the (then already former) East German Academy of Sciences became part 
of DESY. We are glad that this conference attracts since then regulai^ly collegues who contribute 
to the progress in the field. A field, comprising both the branch of applied calculations and that of 
development of new theoretical methods. 

ZFITTER is certainly one of the oldest source-open software projects in elementary particle 
physics with a permanent support. It comprises practically all the theoretical knowledge of rele- 
vance for a precise description of the Z boson resonance in e^e^ annihilation and for Z boson's 

Table 6.2: IBAIKOV=0 (no QCD corrections) or IBAIKOV=2012 [177], AMT4 as described in table 
6.1. The difference to table 6.1: Flag ALEM=2 is chosen with input value Aa^^^{Mz) = 0.02750. 



AMT4 


4 


5 


6 


Diff. 


Exp. Err. 


IBAIKOV=0 




83.9875 


83.9839 


83.9900 


0.0061 


0.086 




2495.2859 


2495.0958 


2495.3662 


0.2704 


2.3 




226.4020 


226.3703 


226.3745 


0.0317 


1.9 


Fw 


2091.1020 


2090.8092 


2090.8474 


0.2928 


42 


Mw 


80.3677 


80.3639 


80.3644 


0.0038 


0.015 


sin^ e'^' 


0.231532 


0.231603 


0.231481 


0.000122 


0.00012 


IBAIKOV=2012 




83.9875 


83.9839 


83.9900 


0.0061 


0.086 


Tz 


2495.0586 


2494.8685 


2495.1389 


0.2704 


2.3 




226.4020 


226.3703 


226.3745 


0.0317 


1.9 


Fw 


2090.8828 


2090.5901 


2090.6283 


0.2927 


42 


Mw 


80.3677 


80.3639 


80.3644 


0.0038 


0.015 


sin^ 


0.231532 


0.231603 


0.231481 


0.000122 


0.00012 



This text is an extended version of the talk. The contribution to the proceedings of LL2012 in "Proceedings of 
Science" (PoS), by A. Akhundov et al., did not appear: PoS(LL2012)036. 
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part in global fits in the Standaid Model [189]. Certainly, today one would create such a project 
quite differently. We can only encourage our collegues to try. Complex projects need (independent) 
duplication. 

Higher order quantum field theoretical predictions face another problem: The solutions be- 
come so lengthy and complex that the idea of source-open software is, in practice, no longer a 
realistic option. This happens already with the ^{ccf) QCD corrections and the complete NNLO 
weak corrections in ZFITTER. They are mere parameterizations of huge, unpublished expressions. 

The LEP/SLC era gave the scientific community unprecedented precision in several fundamen- 
tal quantities like Mz, Tz, the effective weak mixing angle sin^ 6^^, the number of light neutrino 
flavors A^v Of comparable importance is the experimental confirmation of the Standard Model, 
a gauge theory with spontaneous symmetry breaking, as a consistent quantum filed theory, with 
inclusion of higher orders of perturbation theory. 

We are proud that we are being contributing. 
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